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ABSTRACT 

We develop a method for general non-linear cluster lens reconstruction using the 
observable distortion of background galaxies. The distortion measures the combination 
7/(1 — K.) of shear 7 and surface density k. From this we obtain an expression for the 
gradient of log(l — k.) in terms of directly measurable quantities. This allows one to 
reconstruct 1 — /€ up to an arbitrary constant multiplier. Recent work has emphasised 
an ambiguity in the relation between the distortion and 7/(1 — k). Here we show that the 
functional relation depends only on the parity of the images, so if one has data extending 
to large radii, and if the critical lines can be visually identified (as lines along which the 
distortion diverges), this ambiguity is resolved. Moreover, we show that for a generic 
2-dimensional lens it is possible to locally determine the parity from the distortion. The 
arbitrary multiplier, which may in fact take a different value in each region bounded by 
the contour /€ = 1, can be determined by requiring that the mean surface excess vanish 
at large radii and that the gradient of k, should be continuous across k, = 1. We show 
how these ideas might be implemented to reconstruct the surface density, if necessary 
without use of the data in regions where determination of the parity is insecure, and we 
show how one can measure the mass contained within an aperture, again, if necessary, 
without using data within the aperture. 

Subject headings: cosmology: observations - dark matter - gravitational lensing - galaxy 
clusters 
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Observations of giant arcs, arclets and weak dis- 
tortion at large radii show galaxy clusters to be po- 
tent gravitational lenses. With deeper observations 
of very massive clusters it is quite foreseeable that 
we will soon have quite detailed maps of the distor- 
tion field in the best cases, and in both the linear 
and non-linear regime. Here we will address the ques- 
tion: Given only the distortion of background objects, 
how can one reconstruct the surface mass density in 
the lens? The observed surface brightness for an ob- 
ject seen through a lens is /obs(^i) = ftmeii'ijdi)- 
The angular position is measured relative to some ar- 
bitrary fiducial point on the object and where the 
magnification matrix is ipij = 1 — (fi^ij where the 
surface potential (f) is related to the surface density 
by V^(/) = 2S/Scrit = 2k. For a fiat, critical den- 
sity universe Scrit = (47ra;W;(l — wi / Ws))~^ , where 
w = \ — l/i/l + We will assume for simplicity 
that the sources are much more distant than the lens 
{wg ^ w;) so the distortion is effectively independent 
of source distance. This is not really a limitation on 
the method; when this inequality is not strongly sat- 
isfied we actually have more information at our dis- 
posal, but this can be reduced, in principle, to give 
the distortion on a single plane. 

There are many possible statistics which can be 
used to measure the distortion; typically these involve 
measuring quadrupole moments and forming some 2- 
component measure of the polarisation of the image 
shapes. The application of a locally constant shear 
will distort the distribution of observed polarisations 
from its intrinsically isotropic form, so with suffi- 
ciently numerous background galaxies one can mea- 
sure the distortion strength with high precision. The 
best one can hope to learn from such studies is the 
orientation and the ratio of the absolute values of the 
eigenvalues of tpij . This is illustrated in figure 1 which 
shows the distortion of intrinsically circular objects by 
a simple model lens. For concreteness let Lpo be the 
angle of the major axis relative to the a;-axis and let 
i? < 1 be the ratio of the short to the long axis. 

Let us now assume that we are supplied with 
sufficiently precise and detailed maps of R{0) and 
Lpi:j{9). Let the eigenvalues of 4>^ij be ordered such 
that Ai > A2 and let the xi direction lie at an angle 
Lp^ relative to the cartesian a;-axis. In the diagonal 
frame we have tpn = 1 — k — 7 and 1/122 = 1 — k -|- 7 




Fig. 1. — Distortion field for simple "dumbbell" clus- 
ter lens model. At each grid point on the defiector 
plane we have drawn the apparent shape for an in- 
trinsically circular background object. The critical 
lines are shown (solid). The distortion diverges on 
these lines (though note that the arcs are not, in gen- 
eral, aligned parallel to the critical lines). The critical 
density line is shown dashed. The distortion vanishes 
along this line. 
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where 7 = (Ai — A2)/2 is the shear, sc 

7 _ 1 - i)iili)22 



l-K 1 + t/'ll/'/'22 



(1) 



The fact that the observable eigenvalue ratios only de- 
termine this combination of 7 and k was emphasised 
by Schneider and Seitz (1994). They showed that this 
results in a 'global invariance transformation'; there 
is a one parameter family of surface density configura- 
tions which are compatible with any given distortion 
pattern, and we will recover this below. 

Schneider and Seitz also noted an ambiguity in the 
relation between 7/(1 — k) and the observables R, Lpo; 
clearly this depends on the signs of the eigenvalues of 
tpij. Consider first the case where ipiij '4'22 have the 
same sign (even parity image). If they are both pos- 
itive the smallest eigenvalue is ipn, so R = ipiil ''P'n, 
and the long axis is then aligned with the xi direction, 
so Lp^ = Lpcj. If on the other hand they are both neg- 
ative, the smallest absolute eigenvalue is now ip'n, so 
R = lij'nli'w and the long axis is perpendicular to the 
x\ direction so = 1^0 + 7r/2. These appear to give 
different functions for 7/(1 — k) in terms of the observ- 
ables i?, i^o- However, this is in the diagonal frame. In 
the general coordinate frame we have 7]^ = 7 cos 2i^^, 
72 = 7 sin 2i^^ and the numerical value of these com- 
ponents for given i?, 1^0 are the same in both cases, 
so we have 
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even parity 



(2) 



with t\ = ecos2i^o, £2 = esin2i^o and with e = (1 — 
R)/(l + R). The boundaries of the even parity regions 
are the 'outer' critical line where ipn vanishes and 
which lies at k < 1, and the 'inner' critical line where 
i/>22 = and which lies at k > 1. 

Now consider the case that the eigenvalues have 
opposite sign (odd parity image) i.e. ipn < 0. As 
before, the orientation depends on which eigenvalue 
is larger. They are equal and opposite along the line 
K = I where the distortion vanishes. (Note that a 
closed loop on which the distortion vanishes can only 
occur in the odd parity region; in the even parity 
regime the distortion will generally only vanish at the 
singular 'umbilic' points where both eigenvalues are 
equal (Berry and Hannay, 1977).) If k < 1 then R = 
-ipii/'4''22 and ip^ = ipo. Otherwise R = -1P22/1P11 
and Lp^ = Lpo -\- t: /2 and in either case we find 
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odd parity 



(3) 



Thus in either case, the quantity we can hope to mea- 
sure is 7/(1 — k); in the even parity case it is equal to 
the 'ellipticity' parameter e < 1, and in the odd case 
it is equal to the inverse ellipticity. Given perfect ob- 
servations one could simply identify the critical lines 
as lines connecting strongly distorted galaxies and, 
provided the data extend to large radii where it is 
safe to assume the parity is even, one could then un- 
ambiguously determine 7/(1 — k). As we will shortly 
see, for a general two-dimensional lens one can also 
directly determine the parity from local observations 
of the distortion. 

We would now like to eliminate the shear from 
equations 2, 3 to obtain an expression for k in terms of 
observables. Unfortunately there is no unique relation 
between k and 7, but there is a simple relation be- 
tween the gradients of k and of 7: Since ■ji = {{4>,ii — 
22)72, (f)x2] we find d-^i/dxi +d'^2/dx2 = (?i,iii/2 - 
221/2 -I- (/)_i22 = dn/dxi. Similarly, d'^2/dxi - 
d'ji/dx2 = dK/dx2 or 



di 



(4) 



where di 
operator 



d/dxi and where we have defined the 
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d2 
di 



(5) 



It may seem somewhat mysterious that one can 
construct the vector diK out of ■ji (which, while a 2- 
component entity, does not transform as a vector). 
The relation (4) relies on the fact that the shear field 
is derived from a single scalar function 4>, and thus is 
rather special. Replacing ■ji in equation 4 by {1 — K)ei 
we have 



di K — (^1 Ki 'j D ij 6- j 6- j D ij K 



(6) 



or equivalently 

di log(l - k) = -M'^DijCm = Ui (7) 



where 



Mi. 



1 + ei 62 
62 1 - ei 



(8) 



Similar expressions can be obtained for the odd parity 
case (or one can use these expression with e l/e, 
though this will involve differentiating the inverse el- 
lipticity which is probably not such a good idea). 

Equation 7 is the main result of this letter; it gives 
the gradient of the scalar function log(l — k) in terms 
of directly measurable quantities. With sufficiently 
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detailed observations one can simply integrate this 
vector and reconstruct 1 — k, though only up to an un- 
known multiplicative constant (this is Schneider and 
Seitz's global invariance transformation). In fact, the 
multiplicative factor can be set separately in each re- 
gion bounded by the contour k = 1, but this apparent 
ambiguity is removed by requiring continuity of Vk 
across k = I. 

We can now see how one can locally determine the 
parity: the vector u is the gradient of a scalar, so 
curlt/ = Ml 2 — M2,i should vanish, but will typically 
only do so if one is using the appropriate expression 
for 7/(1 — k) in terms of e. This is again because 
the distortion ei(6), while apparently a two compo- 
nent object, has only a single scalar degree of free- 
dom. This introduces special interrelations between 
the distortion values, and this will, in general, not be 
invariant under the transformation e 1/e. This is 
illustrated in figures 2,3 for our simple model lens. 

With the parity determined our problem is for- 
mally solved: one can then determine the difference in 
log(l — k) between any two points by performing a line 
integral. One could, for instance, estimate log(l — k) 
at a point (relative to the very small average value 
around some boundary at large distance) by averag- 
ing over radial lines. In the limit that the shear is 
weak this becomes identical to the method of Kaiser 
and Squires (1993; hereafter KS), and provides a nat- 
ural generalisation of the KS method to the non-linear 
situation. 

It may be that secure identification of the critical 
lines may not be possible. In that case one might still 
be able to use the even parity curl-M statistic as a di- 
agnostic to identify at least roughly the outer critical 
line, and one can simply modify the generalised KS 
statistic to avoid the use of data within the critical 
line (e.g. by averaging only over radial lines which do 
not cross the outer critical line). In this way one can 
establish the surface density k at all points outside the 
outer critical line. This is interesting, but one might 
be more interested in establishing the mass within the 
critical region or within some region with boundary 
outside the critical line. This can be done, and is ex- 
plicitly independent of any data within the boundary, 
as we now show. 

To determine the mass within some aperture we 
need to first determine the shear outside that aper- 
ture. To see this, consider an arbitrary circular 
loop and calculate the mean tangential shear (defined 




Fig. 2. — The ink density is proportional to the abso- 
lute value of curl-M calculated with even parity expres- 
sion for u(R, ipo). The curl vanishes outside the outer 
critical line and in the two regions which lie inside 
the inner critical line, and is non-zero elsewhere, and 
is. Curl-M is particularly large on the critical lines. 
The curl is implemented as a discrete differencing on 
a grid of u values, which themselves involve discrete 
differencing if the gridded e,- values. This causes some 
smearing in the vicinity of the critical lines. 
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Fig. 3. — Curl u calculated with odd parity form for 
u(R,ipo)- The curl is non-zero in the even-parity re- 
gions. It is zero in most of the region lying between 
the critical lines (as it should be), but is also non-zero 
along the critical density contour where u diverges. 
Interestingly, a special case where the this method 
fails is a circular lens, where u is radially directed, 
though the curl-M statistic still responds strongly to 
the presence of the critical lines. 



to be jT = 7i cos -|- 72 sin 2i^) around the loop: 
{jt) = / dip/(2'iT)jT- If we construct local cartesian 
coordinates with radially directed component r and 
azimuthal component /, we have 7^ = ((fi rr—'t'Ji)/'^ = 
(ji^rr — K, SO (jt) = (d/dr) f dip l{2i:)(f) r — (k). The in- 
tegral here is, by Gauss' law simply the mass enclosed 
within the loop divided by the loop circumference and 
we find 



where k is the mean surface density within the loop. 
If we perform the integral J^^ d\nr{jT) (which may 
be expressed as a two-dimensional area integral and 
hence as a sum over discrete data values) we obtain 
the difference between the mean surface density in 
the disk r < ri relative to that in the control annulus 
ri < r < r2 which provides a lower bound on the 
mass interior to ri, tending rapidly to the true mass 
for large r2, while using only data exterior to ri. This 
statistic was used by Fahlman ei al. 1994, and the 
idea is readily generalised to an aperture of arbitrary 
shape. 

How then do we determine the shear in our annu- 
lus? For a sufficiently compact lens k will be negligi- 
bly small in the annulus and we simply use ■ji = e,- 
(and ji = Bj/e^ in the odd parity region). In the 
general case we must first use u to calculate k{9) as 
outlined above (if the annulus were indeed empty we 
would find M = there), and then use ■ji = (1 — K)ei. 

We have described a method which, given suffi- 
ciently precise data, allows one to reconstruct the 
full 2-dimensional surface density for an arbitrary 
lens from measurements of the distortion. The focus 
has been on the mathematics, rather than the 'engi- 
neering' problem of actually measuring the distortion 
field. Our main result is an explicit expression for 
the vector u = Vlog(l — k) in terms of observables, 
and we have shown how this provides a natural non- 
linear generalisation of the KS algorithm. In order 
to recover a detailed picture of the surface density 
in the very centre of super-critical lenses it will be 
necessary to locate the critical lines and apply the 
transformation e 1/e to the data in the odd par- 
ity regime. The odd and even parity curlu statistics 
may prove helpful in this regard as they seem, for our 
idealised data at least, to cleanly detect these lines. 
It remains to be seen whether this will be practical 
with real data. A somewhat less ambitious approach 
is to use the even parity curl-M statistic to identify, 
at least approximately, the odd parity region, and we 
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have shown how the generalised KS method may be 
implemented without using these data and we have 
also shown how the mass within some aperture lying 
outside the critical lines may be measured using only 
the exterior data. 
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